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We establish the set of independent variables suitable to monitor the complicated evolution of the 
spinning compact binary during the inspiral. Our approach is valid up to the second post-Newtonian 
order, including leading order spin-orbit, spin-spin and mass quadrupol-mass monopole effects, 
for generic (noncircular, nonspherical) orbits. Then we analyze the conservative spin dynamics 
in terms of these variables. We prove that the only binary black hole configuration allowing for 
, spin precessions with equal angular velocities about a common intantaneous axis roughly aligned 

to the normal of the osculating orbit, is the equal mass and parallel (aligned or antialigned) spin 
configuration. This analytic result puts limitations on what particular configurations can be selected 
in numerical investigations of compact binary evolutions, even in those including only the last orbits 
of the inspiral. 



I. INTRODUCTION 

^ i' Compact objects are characterized by their size and gravitational radius being comparable. They appear either 
I ' as the end state of the stellar evolution as neutron stars or black holes with a few solar masses (M©) or emerge 
J~ ' ■ from cosmological evolution by continued accretion and a sequence of mergers [l[ as supermassive black holes of 
i ' 3 x 10 6 -j- 3 x 10 9 M , residing in the centers of galaxies. Not much evidence has been gathered for the existence of 
intermediate mass black holes (IMBH), although a detection of a variable X-ray source of over 500 M Q in the galaxy 
ESO 243-49 has been recently reported and interpreted as IMBH 0]. It has been proposed that IMBHs ought to be 
searched for in globular clusters that can be fitted well by medium-concentration King models 

Compact objects are expected to frequently coexist in binary systems, formed either by evolution of a stellar binary, 
by capture events or accompanying the process of galaxy mergers. According to general relativity, compact binaries 
radiate away gravitational waves, a process leading eventually to their merger. Stellar mass binaries are among the 
most prominent sources for the Earth-based gravitational wave detectors LIGO and Virgo [H , while the gravitational 
waves produced during the (low mass) galactic black hole mergers will be sought for by the long-planned space mission 
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The merging process can be split into three distinct phases. By definition the inspiral is the regime of orbital 
evolution, which can be described accurately in terms of a post-Newtonian (PN) expansion. Provided the orbits are 
not excessively eccentric, the same PN parameter characterizes both weak gravity and non-relativistic motion: 
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A manifestly convergent and finite procedure for calculating gravitational radiation to arbitrary orders in a PN 
expansion was proposed [6], based on solving a flat-spacetime wave equation (representing Einstein equations with 
a harmonic gauge condition) as a retarded integral over the past null cone of the chosen field point. A study of the 
Cauchy convergence for PN templates shows an oscillatory behavior: increasing the PN order will not necessarily result 
in a better template Q (2PN templates being closer to numerical results, than their 2.5 counterparts). The predictions 
of various PN approximants (adiabatic Taylor, Pade models, non-adiabatic effective-one-body models) show that their 
convergence to numerical results is comparable It is also known, that alternative template families based on the 
shifted Chebyshev polynomials could exhibit faster Cauchy convergence, than PN templates @. Comparisons with 
full general relativistic numerical runs confirmed that a third PN order approach can be considered accurate for all 
practical purposes. The inspiral is followed by the plunge, where a full general relativistic treatment is necessary, and 
can be handled only numerically; and the ringdown, a process during which all physical characteristics of the newly 
formed compact object are radiated away, except mass, spin and possibly electric charge. 

In this paper we investigate the conservative dynamics during the inspiral of a spinning compact binary system. 
We include spin-orbit (SO), spin-spin (SS) and mass quadrupole - mass monopole (QM) couplings, each to leading 
order. The precession due to these interactions was first discussed in [l(|-[llj]. With the spins and mass quadrupole 
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moments included, the number of variables in the configuration space increases considerably, therefore we propose to 
find a minimal and conveniently chosen set of independent variables. 

We note discussions of various aspects of the dynamics and gravitational radiation related to the SO coupling in 
12 |-[l4j], SS coupling in [li{-[l6j|, and QM coupling in (TtJ-QjI- PN corrections to the SO coupling were presented in 
20| and the Hamiltonian approach including spins has been also worked out [21] . Most recently, the back-reaction 
on the dynamics due to asymmetric gravitational wave emission in the spinning case, possibly leading to strong recoil 
effects, has been widely investigated, both analytically [l4j |. [22j and numerically for particular spin configurations 
[23| . Empirical formulae giving the "final spin" have been advanced in Refs. [2~ij and some of them compared in [25| . 
Zoom- whirl orbits (generic for particles orbiting Kerr black holes (26[) were also found in the framework of the PN 
formalism 27]. A larger spin increases the likeliness of apparition of such orbits (28|. Gravitational wave emission is 
hold responsible for the occurrence of the spin-flip phenomenon [29]- 30] in X-shaped radio galaxies [29], 31]. Recently 
it has been shown, that for a typical merger of mass ratio at about 0.1 the combined effect of SO precession and 
gravitational radiation will result in the spin- flip occurring during the inspiral [32j . 

In Sec. [n] we introduce the set of dynamical and configurational variables characterizing the compact spinning 
binary. Both the configurational and a subset of the dynamical variables depend on the choice of the reference 
system. We use a number of four such systems, to be defined in subsection III Bt only one of them inertial, the rest 
of three being rather adapted to the binary configuration. In subsection III CI we derive two relations among the time 
derivatives of the introduced angular variables. As a result the time evolution of the configurational variables is 
determined by the evolution of one single configurational angle a and the true anomaly Xv At the end of this section 
we express the position and velocity vectors in the chosen reference systems. As a by-product we recover the true 
anomaly parametrization of the radial evolution, valid for the chosen perturbed Keplerian setup. 

Sec. Mil introduces the angles characterizing the angular momenta (total and orbital angular momenta and spins). 
The number of independent variables characterizing them is shown to be 6. We will chose them either as 5 angles 
and a scale, or equivalently as 3 angles and 3 scales. 

In Sec. IIVI we analyze the conservative evolution of the spins, which is purely precessional, with the inclusion of 
the leading order spin-orbit, spin-spin and mass quadrupole - mass monopole couplings. We clarify the order (both 
PN and in the mass ratio) at which the various contributions occur. Then we investigate, whether there are spin 
configurations conserved by precessions, and we derive a no-go result. 

The gravitational constant G and speed of light c are kept in all expressions. For any vector V we denote its 
Euclidean magnitude by V and its direction by V. 

II. KINEMATIC AL AND DYNAMICAL VARIABLES 
A. Variables 

We consider three distinct set of variables. 

(a) The physical parameters of the binary: The two compact objects are characterized by masses m„ spins Sj 
(i = 1,2), and mass quadrupole moments. 

Equivalently to m, we can use the total mass m = mi -I- TO2 and the reduced mass /1 = raxmijm. We assume that 
mi > 77i2. We also introduce the mass ratio v = m^/mi < 1 and the symmetric mass ratio 77 = \xjvn 6 [0,0.25]. The 
two mass ratios are related as 

v , > 

V = 2 . 2 

and for small v we have r\ = v — 2v 2 + O (^ 3 ) . We also note the useful relations 

2 2 2i— 3 / \ 

m i =mr\v . (3) 

Equivalently to Si we can introduce their magnitude, polar and azimuthal angles. It is convenient to define 
dimensionless spin magnitudes \i G [0, 1] by 

Si = —m 2 lXl = —m 2 r)v 2i - 3 Xi • (4) 
c c 

As for the spin angles, they depend on the chosen reference system. We will discuss various possibilities in detail in 
Section Unl 
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We consider axisymmetric compact objects. Therefore the mass quadrupole of the i axially symmetric binary 
component is characterized by a single quantity Qi, its quadrupole-moment scalar [TtJ- Provided the quadrupole 
moment originates entirely in its rotation (what we shall assume), then the symmetry axis is and 

G 2 

Qi = T w ^ m \ > ( 5 ) 

with the parameter w S (4, 8) for neutron stars, dependin g on their equation of state, stiffer equations of state giving 
larger values of w [33| , [TtJ ■ For rotating black holes w — 1 |34j . The negative sign arises because the rotating compact 
object is centrifugally flattened, becoming an oblate spheroid. 

(b) Dynamical variables: Up to 2PN accuracy the energy E and the total angular momentum vector J = L + S1 + S2 
are conserved [12j . The orbital angular momentum L and the spins Si are not conserved separately, as the spins 
undergo a precessional motion. This will be discussed in detail in Section IIVI 

(c) Angular variables characterizing the orbit: 

The instantaneous orbital plane is perpendicular by definition to the Newtonian orbital angular momentum Ln = 
jir x v and it evolves due to the spin precessions. We define (i) the inclination a of the orbital plane with respect 
to the plane perpendicular to J (thus a is the angle span by Ln and J); (ii) the angle 4> n between the intersection 1 
of these two planes and an (arbitrary) inertial x-axis x taken in the plane perpendicular to J, finally (iii) the angle 
tjjp measured from 1 to the periastron (see Figs Q] and [3J the indices p and n stand for the periastron and node line, 
respectively). 




FIG. 1: The polar and azimuthal angles of the total angular momentum JJ, Newtonian orbital angular momentum LjvLn and 
spins Si, 281,2- Azimuthal angles are shown in the non-inertial system ICj = ^1, k, J j , polar angles both in /Cj and respective 

to Ln- The relative angle of the spins is 7. The non-inertial character of the system /Cj is encoded in the evolution of the 
angle cj> n , measuring the angular separation of an inertial axis x from the axis 1. 



B. Reference systems 



For a better bookkeeping we introduce the inertial system /Q with x and J standing as the x- and z-axes and three 
non-inertial systems ICj, Kl and K-a- 
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FIG. 2: The relative angles of the total angular momentum JJ, Newtonian orbital angular momentum LjvLn and spins 
Si,2Si,a as in Fig[T] The intersection of the planes perpendicular to Ln and J, respectively is the node line 1. The inertial 
axis x is at angle (j>„, measured from 1 in the plane perpendicular to J. The azimuthal angles (ipi, tp2, ip P ) of the spins and 
Newtonian Laplace- Runge-Lenz vector tIjvAn (pointing towards the periastron) are also measured from 1, however in the 
plane perpendicular to Ln. The true anomaly \v i s t ne angle between An and the position vector rf. Two of the basis 
vectors of the inertial reference system /Q = ^x, y, J^j and of each of the three non-inertial reference systems /Cj = ^1, k, J | 

Kt 



(l, rh, L N ^ , K-a = (A N , Q N , L N j are shown. 



In the system KLj the 2-axis is fixed along J, while in /Cz, along Ln- We choose l = Jx LN/sina as the x-axis of 
both systems. The system ICj is complete by k = J x 1 and JCl by m = Ln X I- 

The system Ka also has Ln as the z-axis, however its x-axis is defined by the Laplace-Runge-Lenz vector 



An = v x L 



N 



Gmfj, 



(6) 



which satisfies the constraints 



A 2 



2EnL 2 v 
A 4 



(7) 



and Ln-An = 0. Here r and v are the position vector and velocity of the reduced mass particle orbiting m. The 
y-axis is defined by Qn = Ln x An- The orthonormal basis of K-a is therefore (An, Qn, Ln). 

The three angles ((p n ,a,ip p ) will be referred to occasionally as Euler angles, as three consecutive rotations with 
—(f> n , & and ip p about the axes z, x and again z transform as /Q — > JCj — > Kl — > JCa- The sequence of these rotations 
is encompassed in the transformation matrix 



U(-(p n ,a,ip p ) = U z (ip p )1Z x (a)TZ z (- 

scj) n + sin tp p 



cos ipp cos <p n + sin tp p cos a sin <p n — cos tp p sin 
- sin ip p cos <pn + cos ip p cos a sin <p n sin tp p sin <p 
— sin a sin <p n 



t> n + sin ipp cos a cos <p n sin tp p sin a \ 

- cos ipp cos a cos <p n cos tpp sin a , (8) 
cos a I 



— sin a cos (p n 



where 1Z with one argument denotes the corresponding rotation matrices acting on the coordinates. 
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C. Constraints on the Euler angle evolutions 

The coordinates of the reduced mass particle in the inertial system /Q can be obtained by applying the transfor- 
mation 1Z (— 0, —a, <fi n ) to the coordinates of the vector r = r (1, 0, 0). Here = ip p + \p is the angle span by 1 and 
f , with Xp defined as the true anomaly, the angle span by An and f . We obtain 

cos <fi n cos + sin <fi n cos a sin 

— sin (fi n cos + cos <fi n cos a sin | . (9) 
sin a. sin ip 

A tedious, but straightforward computation carried on in the system Ki gives 

k / sin a sin [cos (fi n cos ip + sin <fi n cos a sin 0] 

— — = (fi n I sin a sin -0 [— sin <fi n cos -0 + cos (fi n cos a sin 0] 
//r \ sin 2 a sin 2 0-1 

— sin 0„ cos a cos + cos <fi n sin 
+d sin ip | — cos <fi n cos a cos — sin <fi n sin 
— sin a cos 

Sin (fin sin a 

+0 I — cos</>„sina I . (10) 



cos a 



From here we readily obtain 



T 2 

2 4 

fi z r^ 



^0 — n cos cxj + (j> n sin a cos + dsin0^ . (11) 



Also, dividing the third component (which by definition is (L\) z = -L^cosa) by cos a we get 

— ^ = if b — (fin cos a — (^fi n sin a cos + d sin 0^ tan a cos . (12) 

In the Newtonian approximation the Euler angles being constant, we recover — fir 2 x P and (Ln) z = fir 2 x P cos a. 
By squaring Eq. (fT2"|) and subtracting from Eq. (jTTJ) we obtain the identity: 







fl - tan 2 



a cos 2 0) ^</>„ sin a cos + d sin 0^+2 tan a cos ^0 — (fi n cos 
x (^<fi n sin a cos + d sin 0J . (13) 

The first factor cannot vanish, as to Newtonian order it gives 2 tan a cos tfiLff/fir 2 ^ 0, therefore the vanishing of the 
second factor (by reintroducing = p + \p) gives: 



• . tan (0 p + x P ) ,,„s 

0„ = -a : . (14) 

sma 

Reinserting this in either of the Eqs. (fTl) or (|T2l gives 

P + X P = — % + (fincosa . (15) 
fir 2 - 

We have just derived two relations among the time derivatives of the Euler angles and of the true anomaly, which 
restrict the number of independent angular variables introduced up to now to a and Xp- 

D. The position and velocity vectors in the bases K,a and ICl 



Simple computation starting from the definitions of An and Qn gives 

(2E N Gm\ 

An = M 1 r—firrv 

V fi r J 



Qn = Gmfi 2 rr+ (L 2 N - Gnifi 2 r) v . (16) 



From here the expressions of the position and velocity vectors in ICa emerge as 



r = A N + — rrQ N , 

Gm/i . -j L N ( 2E N Gm\ A 
v = t — rA N + — + Qn • 

A N A N \ n r J 

In terms of the true anomaly Xp (the azimuthal angle of r in the system ICa), the position vector is given by 

r = r (cos x P A N + sin x p Qn , 

which compared with Eq. (|17p gives the true anomaly parametrization: 

T 2 



/i (Gm/i + Aat cos Xp) 
A N . 

— smxp ■ 

Ln 



In terms of the true anomaly, the velocity is expressed as 

Gm/i 



L 



N 



sinxpAN + cosxp 



A 



N 



Grrifi 



Qn 



Its square gives v 2 in terms of the true anomaly: 



(Gm/i) 2 + A 2 N + 2GmfiAN cosx P 



T 2 



(17) 
(18) 

(19) 

(20) 
(21) 

(22) 
(23) 



(The same emerges from the definition of the Newtonian energy En = /iu 2 /2 — Gm^i/r, by applying Eqs. ([7]) and 
As the basis vectors of ICa are related to the basis vectors of ICl by a rotation with angle ip p : 



An = cos ip p l + sin -0pm , 
Qn = — sin^pl + cos V>p m 



it is straightforward to rewrite r and v in the basis ICl as 



r = r (^cos ijA + sin ijjrh ) , 
1 



(Gm/i sin if> + An sinV>p) 1 + (Gm/i cos ^ + An cosip p ) m 



(24) 

(25) 
(26) 



III. CONSTRAINTS ON ANGULAR MOMENTUM VARIABLES 
A. The 5 angular degrees of freedom 

The polar and azimuthal angles of Ln and Si in ICj are (a, — tt/2) and (/3i, </>i), respectively, such that 

Ln = — sin ak + cos a J , 
Si = sin/3j cos 0il + sin/?i sin + cos/?^ J . 

Similarly, the polar and azimuthal angles of J and Si in ICl are (a, tt/2) and (Kj,^>i), respectively, thus 

J = sin am + cos oLn , 

Si = sin Ki cos tpi\ + sin m snn/>irii-|- cos k^Ln . 



(27) 
(28) 



(29) 
(30) 
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By comparing the two forms of the 1 component of the vectors Sj we get 

sin Ki cos ipi = sin /3j cos (pi . (31) 

By computing Si • L N in both systems we find 

cos Ki= cos a cos /3j — sin a sin ft sin <pi . (32) 

As shift = — cos (7r/2 + ft) and 7r/2 + ft is the relative azimuthal angle of Ln and Si, Eq. (|32|) is but the spherical 
cosine identity in the triangle defined by these three vectors on the unit sphere. 

Similarly, from the two expressions Si - S2 = cos 7 written in both reference systems we find the spherical cosine 
identities: 

cos 7 = coski COSK2 + sinKi sin«2 cos Aip , (33) 
cos 7 = cos ft cos ft + sin ft sin ft cos Atp . (34) 

where Aip = ip2 — V'l and A<p = ft — ft are the differences in the azimuthal angles of the two spins in the two systems 
JCl and JCj, respectively. 

Other spherical triangle identities arise by computing Si • J in both systems: 

cos ft = cos a cos Ki + sin a sin Ki sin ipi . (35) 

Then Eqs. pip and (|35|) give Pi, <pi as function of Ki,ipi and a. Inserting these in Eqs. (1331) and (|34p and eliminating 
7 could in principle give a as function of Ki,ipi alone. We get: 

sin /3i sin <pi sin P2 sin ^2 = (sin a cos Ki — cos a sin ki sin ) (sin a cos K2 — cos a sin K2 sin ip2 ) ■ (36) 

As the orientation of the spins are independent, we obtain 

sin Pi sin <pi = sin a cos Ki — cos a sin sin ipi , (37) 

however the direct computation of the left hand side by employing Eqs. (f3"T|) and (|35|) results in the right hand side, 
leading to an identity rather than an expression of a as function of Ki,ipi. Therefore Eq. (|34j) can be considered as a 
consequence of the other equations. Similarly one can show that Eqs. (|32[) are consequences of the other equations. 

We conclude that there are 5 independent constraint equations for the 10 angles (a, ft, </>,-, k,, ipi, 7), namely Eqs. 
(|3"Tj) . (|35|) and (|33[) . and we can take (a,Ki,tpi) as the independent angles. The network of all angles in the systems 
Kj and Kl is represented on Figs [T] and [SJ respectively. 



B. Orbital angular momentum 

The total orbital angular momentum L contains pure general relativistic (PN, 2PN) and spin-orbit (SO) contribu- 
tions [HI: 1 

L = Ln + Lpn + Lso + L2PN • (38) 

There are no spin-spin or quadrupole-monopole contributions to the orbital angular momentum [l8j . Here the Lpn 
and L2PN contributions are aligned to L N (cf. Eq. (2.9) of Ref. (H): 

Lpn = epatLn , 

1 — 377 /v\ 2 . Gm 



- + 3 + 7/ — , 39 



The equations of motion leading to this expression were derived in harmonic coordinates, imposing the covariant spin supplementary 
condition. 
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and 



1*2PN — £2P7vLn , 

E2ra = 5(i-T, + i3fl(H) 4 -I,(2 + 5,)f=(r 

4p-^-V)^© 2 + i(i4-4i, + V)(^)' 

The SO contribution (Eq. (2.9.c) of Ref. [lj]) can be rewritten as 

2 

Lso - S > [ £ i ( f ' ^i) ? + (v • Si) v- (e? + ej) S ; 

Gra 



(40) 



2 

,3-2 



I- 2 



Note that 



and 



e PW = 0(e) 0(1, r?) , 

£2PA , = O^Oil^.if) , (42) 



e[ = 0(e)0( n )0(l^) , 

(? = O (e) O (rj) O (v 3 - 21 ) . (43) 

In order to evaluate the PN order of the Lso contribution in J, we evaluate on circular orbits 

Si_ = (G/c)m 2 W 2 ^ Xl = /Gm\ £^- 3 
Ljv firv \ c 2 r J v 

= 0(j' a )0 Xi , (44) 

which continue to approximately hold for eccentric orbits. This reasoning shows that the SO contribution is of 1.5 
PN order and also indicates how to pick up the dominant terms when the mass ratio is small or when one would like 
to employ a less accurate, but simpler description, dropping higher order terms. 
The total angular momentum is then 

2 

JJ = (1 + e PN + e 2PN ) L N t N +Y, S i [ £ i ( f ' Si) f + e? (v ■ S ; ) v+ (1-ef - eft Si] . (45) 

i=l 

C. One scaling degree of freedom 

In this subsection we will employ the projections of the Eq. (|45|) in order to derive relations between the angles 
and magnitudes of the angular momenta involved. In the ICl system the projections along the axes 1, m and L\ give, 
respectively: 

2 

= Si sin m {(1 — ef — e-') cos + e[ cos (ip — ipj) cos?/' + elfS [GrnfJ, sin (ip — ipj) + A N sin (ip p — ipi)]} (46) 

i=l 
2 

Jsina = ^ Si sinKj {(1— e[ — £j ) sin^ + e\ cos (ip — ipj) simp — eJC [Gm^sin (ip — tpj) + sin (ip p — ipj))} (47) 

i=l 

2 

J cos a = L N (1 + e PN + e 2 pjv) + 2J Si (1-ef - eJO cos K * j ( 48 ) 
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where 



Grrifi sin (ip p + Xp) + An sin tp p 



Gm[i cos (ip p ^ 



- 2Gm^Ajq cos x P 
X P ) + A N cos ipp 



In the derivation of Eqs. (146 



(Gm/i) 

48l) we have employed Eqs 



2Gra[iAN cos Xp 

(l26l). (l29l). (l30l from where we also obtained 



f • Si = sin Ki cos (V> 

~ sin k 
v • a, = 



L N v 



—Gra[i sin (i/> — tpi) + An sin — ip p )] , 



(49) 



(50) 
(51) 



with v given by Eq. (|23|) . 

We thus have introduced the 14 quantities ( J, L, Xi, A; <Pi> K ii ifri> l) describing the angular momenta, which are 
constrained by 8 independent relations. This leaves us with 6 independent variables. 5 of these can be thought as 
the angles defining the directions of the spins and orbital angular momentum in the Kl system [a, Ki,ipi), a sixth 
one being a linear scale, most conveniently chosen as J. 

Note that in Eqs. (f46|) - (|48|) the coefficients epn, £2PN, ef, ^ depend only on the masses and Xv Therefore all 
dependences on ipi are explicit. In principle Eqs. (|46l) - (jT7|) can be used to express tpi as function of Ki, a, ip p , the 
masses and the spins Xi- ln practice however this may be cumbersome. The easiest way to do it is to rewrite both the 
sin^i and cosf/'i in terms of the variables Xi — tan^i/2. Then Eqs. ([4^1) - (j4T)) become second rank coupled polynomial 
equations, possibly leading to two distinct values of tpi for each Xi- 

Finally, Eq. (|48p can be employed to eliminate Ln in the detriment of the angular variables, spins and J, by a 
series expansion in e to 2PN order accuracy as 

2 

Ln = J (l - epw + £ 2 PN - £2Pn) cos a - ^ S. t (l - e PN - e\ - e? + e PN (\ + e PN e\ + e 2 PN - e 2PN ) cos 

i=l 

2 

= (l - e PN + £ 2 PN - £2Pn) L N: o + ^S t (I - e PN ) {t\ + e-i) cos Ki , (52) 



i=i 



where 



Ln,o — J cos a — Si cos Ki 



(53) 



is the leading order contribution to the orbital angular momentum, arising when we approximate J as the sum of the 
Newtonian orbital angular momentum and the spins. For convenience we also give 



1 _ 1 + £pn + £2Pjv J2i=i & (1 + epjy) (ef + eX) cos Kj 
Ln 



L 



JV,0 



T 2 



(54) 



D. Summary: the independent variables 

The considerations in this section leave us with the following alternative sets of independent variables, all charac- 
terizing the angular momenta: (a, Ki,tpi, J) or (a, «;,-, J). The second set represents the most advantageous way of 
choosing the variables. Most notably, while tpi are constant over the orbital scale, they vary with the precessions. By 
contrast Xi are constant over the precession time-scale either, moreover they are unaffected by gravitational radiation 
reaction, to quite high PN orders. Also, J stays constant up to 2PN accuracy (thus over precession time-scale) as 
opposed to either of L, Ln, Ln,o- It changes only over the radiation time-scale. 

Once the evolution of Xp is known, the other two Euler angles [4>n,'4 , p) become determined by the rest of variables 
through Eqs. (H]) and (fl%|) . 
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IV. SPIN EVOLUTION 



The spins obey a precessional motion, as was derived for bodies with arbitrary, but constant mass, spin and 
quadrupole moments (Eqs. (39) and (43) of Ref. [Rj, see also Ref. (llj): 



(55) 



with the angular velocities consisting of SO, SS and QM contributions. The latter come from regarding each of the 
binary components as a mass monopole moving in the quadrupolar field of the other component. 
The precessional angular velocity is decomposed as 



S7i 

nf° 
n9 M 



GU + Zv 



3-2i\ 



2c 2 r 3 



-L Ln 

N 



3(f ■ Sj) f - S 



3Gm,jQi 
r 3 S, 



r -S; f 



where j ^ i. The sum of the SS and QM contributions, by employing Eqs. (|3|)-([5]) is 

G 



3(f -Sj)f + 3w 2(i ~ i) (f SOf-Sj , j^i. 



(56) 



(57) 



In order to evaluate the PN order of the coefficients in Eqs. (|56l) . we will employ the estimate from a footnote of 
Ref. [R|, according to which 



0{T- r ) 



(58) 



T being the radial period, defined as twice the time elapsed between consecutive f = configurations. We obtain 



G (4 + 3z/ 



3-2i 



1 Gto v c L 



N 



2c 2 r 3 



N 

GSi 



2 c 2 r c r mrw 



4 + 3z/ 



0(e)0(v- 1 ,l,v)0{T- 



3GrrijQj ( Gm\ c „ / o/ \ „ , v „,„ 
3^ = 3w — -r?Xi = U 3/2 O WxiO r 1 



(59) 



Thus on the orbital time-scale the SO precession is a 1PN effect, while the SS and QM contributions appear as 1.5 
PN corrections. As both the SO and SS angular velocities contain terms with O {y^ 1 ) factors, whenever the mass 
ratio is small, the respective precessions amplify. 

As ttf M oc Xii the QM precession qualifies as a self-spin effect. 



A. Spin configurations preserved by precessions 

With only the leading order SO precession taken into account, both spin vectors undergo a precession about Ln- If 
777-2 = also holds, the instantaneous angular velocities of the precessions are identical, and the spin configuration 
is preserved with respect to the osculating plane of the orbit, rigidly rotating about its normal. 

With the SS and QM contributions to the spin dynamics included, the above simple picture does not hold any 
more. In the remaining part of this section we analyze whether there are spin configurations which are preserved by 
precessions, in the sense that they rigidly precess about a common rotation axis. 

We will carry on this analysis order by order, starting with the leading order SO precession. One possibility is that 
both spins are either aligned or antialigned with the orbital angular momentum Sj = ±Ln, then there is no precession 
at SO order. Moreover, at the next order we immediately obtain fl9 M = and ttf s oc Ln, such that Si = 0. Thus, 
when the spins are perpendicular to the osculating orbit at some initial instant, they stay so, even with the SS and 
QM parts of the dynamics included. 
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Another possibility to consider is, that the two spins precess with the same angular velocity about a common 
axis. We could check, whether the axis defined by Eq. (|5B|) could be this, however we will allow for more generic 
possibilities. As Si undergo pure precessions, one can add arbitrary contributions (G/c 2 r 3 ) (Vi — 1) Si to S~2i without 
changing the dynamics, and ask the question, whether a common instantaneous axis of precession exists for both spin 
vectors, about which they precess with equal angular velocities, such that fl^ = fi 2 ? The condition for this would be 



= 



-L N + {f • [(!-■ 



S 2 - (1 - wv) Si] } f-i (P 2 S 2 - PxSO . 



(60) 



For Vi of order unity (meaning that this axis is not very far from the normal to the osculating orbit) the leading order 
contribution in Eq. (|60[) remains the term proportional to Ln, the vanishing of which implies rri2 = ni\. For the next 
order then we get 



P 2 S 2 - 7?! Si = 3 (1 - w) [f • (S 2 - Si)] f 



(61) 



For black holes (w = 1) this gives P2S2 = Pi Si, thus the spins should be parallel (aligned or antialigned), and the 
common axis of synchronous rotation is 



G 



L N + 3(r-S)f-S 



V 1 S 1 +V 2 S^ 



(62) 



with S = Si + S 2 . Neither the axis of rotation nor the angular velocity are unambiguous, as fi[ depend on Vi, 
however the axis stays close to Ln (no choice of Vi would render the axis of rotation exactly to Ln)- In summary, 
only parallel black hole spins can rotate with the same angular velocity about a common axis, provided the axis is 
only slightly different from the normal to the osculating orbit. 



V. CONCLUDING REMARKS 



In this paper we have derived the set of independent variables suitable to monitor the evolution of a compact 
spinning binary during the inspiral. The number of independent variables characterizing the spins and orbital angular 
momentum was shown to be 6. We have chosen them either as 5 angles and a scale, or alternatively as 3 angles and 
3 scales. For the first choice we found advantageous to employ the magnitude J of the total angular momentum; 
the angles a and ft, span by the Newtonian orbital angular momentum Ln with the total angular momentum J 
and with the spins, respectively; finally the azimuthal angles ipi of the spins in the plane of motion (perpendicular 
to Ln), measured from the ascending part of the node line (the intersection of the planes perpendicular to J and 
Ln-) For the second choice we propose J, a, n t and the normalized magnitudes of the spins Xi- As both J and 
Xi are unaffected by precessions; moreover Xi vary extremely slowly with gravitational radiation reaction, the latter 
set seems more advantageous. Nevertheless, expressing tpj in the detriment of Xi is not immediate (the respective 
equations are provided). 

These 6 variables have to be supplemented by the true anomaly Xp- The non-inertial character of the reference 
systems introduced in Section 3 can be specified through one single angle <j) n , characterizing the node line, the 
evolution of which is governed by the spin-orbit coupling. The orbital evolution being quasi-Keplerian, the position 
of the periastron is specified by an evolving angle ip p . As shown in subsection III CI the evolution of these two angles 
follow from the evolution of a and Xv 

In this paper we have also proven a no-go result, according to which in a 2PN accurate dynamics, with the leading 
order SO, SS and QM precessions included the only binary black hole configuration allowing for spin precessions with 
equal angular velocities about a common instantaneous axis roughly aligned to the normal to the osculating orbit, 
is the equal mass and parallel (aligned or antialigned) spin configuration. When including only the SO precessions, 
the equality of masses is required, but there is no constraint on the spin orientations. By approaching the innermost 
stable orbit, the PN parameter increases (leading eventually to the breakdown of the PN expansion), such that the 
importance of higher order contributions is enhanced. Therefore the SS and QM precessions (of higher order than the 
SO precession), which lead to the above constraint on the spin directions, become increasingly larger. The result thus 
will hold up to the very last orbits of the inspiral, and to the extent the PN result approximates well dynamics there, 
during the plunge. This analytic result puts limitations on what particular precessing configurations can be selected 
in numerical investigations of compact binary evolutions, even in those including only the last orbits of the inspiral. 
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